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Abstract 


^This irr-n-nlTili-rn Hnrrihrt jtinHrling and controHTesifen for flexible manipulator!, both 
from an experimental and analytical viewpoint/r rom the aplication perspective, 

£ n ongoing effort within the laboratory environment at The Ohio State University, 
where experimentation on a single link flexible arm is underwayr'Several unique features 
of this stmly are described here. First, the manipulator arm is slewed by a direct drive 
dr motor and has a rigid counterbalance appendage. Current experimentation is from two 
viewpoints: 1) rigid body slewing and vibration control via actuation with the hub motor, 
and 21 vibration suppression through the use of structure-mounted proof-mass actuation 
at the tip. Such an application to manipulator control is of interest particularly in design 
of space-based telerobotic control systems, but has received little attention to date. From 
an analytical viewpoint, wgadtafiiua parameter estimation techniques within the closed- 
loop for self tuning adaptive control approaches c"At?(Vititrodured is a control approach 
based on output feedback: and frequency weighting to counteract effects of spillover in 
reduced -order model design. A model of the flexible manipulator baaed on experimental 
measurements is evaluated for such estimation and control approaches. 
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1. Introduction 

Traditionally, robotic manipulator arms have been modeled as being composed of rigid links, with co-located actuators and sensors, 
towards the goal of ensuring stable and reliable control. In order for typical manipulator arms to maintain this rigid property as 
modeled while carrying payloads, the mechanical design requires large and massive links. This in turn dictates that the torques 
applied by the joint actuators be large, and heavy, usually geared motors are needed for actuation. Moreover, the controller for 
such a system is forced to move the arm slowly and deliberately so as to prevent any swaying or vibrations. 

In recent years there has been muck interest in using light-weight, higher performance arms for both commercial and space-baaed 
applications, leading to the study of flexible manipulator control. The advantages of flexible robotic manipulators are many, 
including faster system response anil lower energy consumption, smaller actuators and overall trimmer mechanical design, reduced 
nonlinearity effects due to elimination of gearing, less overall mass and generally less cost. Obvious tradeoffs, however, complicate 
the issue of flexible manipulator control, primarily centering on the design of controllers to compensate for. or to be robust in 
the presence of flexure effects. With the advent of advanced computational resources, strides are currently being made towards 
solution of (he many problems associated with control design. 

Control design for lightweight flexible manipulator arms has gained the attention of control theorists only recently, and several 
approaches have emerged. Most prominant are approaches which either linearize and truncate for controller design, or solve the 
nonlinear robotics problem for rigid link motion control and treat the flexible dynamics separately. For example, the problem of 
observation spillover and truncation error effects is treated in lj, where in simulation studies a linear feedback scheme around 
a reduced order model is introduced for a single-link manipulator. In 2! control of the rigid motion is accomplished via state 
feedback linearization whereas vibrational dynamics are treatrd as disturbance effects. Several other analyses have appeared along 
these hasi. lines, using various approaches ■I.C.Vti. O. From an applications viewpoint, however, only a few studies have been 
documented lor parameter estimation, system identification, and control. Most prominant among these are tbr works of Book. , 1 
<il. H.q.lll. 1 1 for time optimal slew experiments, related studies at JPL in flexible beam control 12.131, Schmitz and Canon I t 
using non colocated and tip position sensing in the control algorithm, and several studies conducted at NASA LaRC 1'i.lfil. 

In this presentation we report on progress made to date on modeling and control design for flexible manipulators, both from an 
experimental and analytical viewpoint. Specifically, we discuss the ongoing effort within the Control Research Laboratory at The 
Ohio State University, where experimentation on a single link flexible arm is underway. The manipulator arm is slewed by a direct 
drive dr motor and lias a rigid counterbalance appendage. Current experimentation is from two viewpoints: 1 1 rigid body slewing 
and vibration control via actuation with the hub motor, and 2) vibration suppression through the use of structure-mounted proof- 
mass actuation at tin- tip. Real-time parameter estimation techniques, within the closed-loop for self-tuning adaptive control, is 
under investigation and is described briefly here. In these initial studies, a model of the flexible manipulator based on experimental 
measurements is evaluated. 
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7 . Experiment*! Setup 
i.l Apparatus 

Within the Control Research Laboratory at The Ohio Stela University, several experimental configuration* art under etudy 
for *y*tem identification and dewing and vibration control for flexibte mechanical itructore*. In thi* presentation we local on 
experimentation and eimulation analyti* with a tingle link flexible arm, depicted in Figure 1. The arm i* made of 0.0925 inch 



Figure 1: OSU Flexible Arm 

thick aluminum and n counterbalanced with a rigid aluminum appendage with mat* equal to that of the arm. Hob actuation i* 
accomplished by a .1. 1 ft - lb direct drive dc editor which ha* an optical encoder with a quadrature digital output to icnie motor shaft 
position. and a tachometer to measure motor shaft speed. This, then, allows both hub position and velocity feedback for cmstrol. 
Strain gauges (for monitoring and parameter estimation) and an accelerometer are placed along the arm. and a 512-elemen» CIO 
linear array camera with RS-422 interface is used for sensing the tip position by observing the lamp fixed to the tip of arm. With 
such a scheme, the tip tensing mechanism (camera) it utilised in verification and tuning of the predicted endpoint poaitioss. A 
related objective for this setup it to achieve control without camera information feedback, with for example rate and acceleration 
sensing feedback, for application in space-based manipulator systems where off-structure reference for sensing is impractical. 

Some characteristics of the arm are given below, 

Thble 1; Arm Characteristics 

! Material 6061 -T6 Aluminum 

l 

| Modulous of Elasticity 68.944 * 10* N/m l 

Cross Sectional Area Moment of Inertia 

Flexible Arm 3.350 < 10*" m 4 

Rigid Appendage 2.427 < 10'* m 4 

Lengths 

Flexible Arm 1.0 m 

Rigid Appendage 0.381 m 

The unique features of the structure are the direct drive mechanism, chosen to minimize effects of backlash and other nonlinearities 
due to gearing, and the counterbalance appendage, which provides a more realistic model of application-oriented structures. Such 
a hybrid structure does, on the other hand, pose unique problem* for analytical modeling. 

Two computing environments are available in the laboratory for real-time control and data acquisition. The first such system is *n 
IBM AT which uses different combinations of several custom-built cards in addition to the A/D equipment. These cards include a 
controller for the slewing motor with electronics utilized in procezzing data received from the linear-array camera. Another card, 
designed in-house, processes strain gauge and accelerometer data, and includes a low-noise, higb-gain amplifier with a low paai 
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filter and excitation to the bridge circuitry. A third cuitom board it used to drive the proof-mazt actuators (discussed in [17| and 
later in this presentation). The board receives an analog voltage from the D/A and amplifies it to a current which is adequate to 
drive the actuators. The linear array camera is interfaced to the computer using a custom board which converts the camera's serial 
data stream to a number corresponding to the position of the beam endpoint. The second computing system is the MicroVax II, 
equipped with commercially available A/D-D/ A boards and real-time operating system software. For the study presented here, 
the data acquisition is carried out using the IBM AT due to the availability of the camera interface electronics. 

!.! Modtltnf and Frtqucncf Reipontt 

For purposes of finite element modeling studies, the arm is enalysed in two separate components, those being the flexible arm 
itself and the rigid counter balance. The cylindrical mass at the end of the rigid component is modeled as a point mass, and 
the two components are connected at the pinned joint (motor shaft). For the FEM analysis, each component is modeled as a 
two-dimensional elastic beam, and the software package ANSYS [18| was used to generate the first five modes of the system 
shown in Table 2 below. We note that torsional modes we re assumed to be insignificant, and were therefore neglected in the 
analysis. The principle advantage to modeling the arm in this manner is that the effects of the counter balance in the static 
characteristics are included. Several other approaches were utilised, such as considering the joint at the motor shaft to be a fixed 
point (clamped freel. negating any effect the counter balance may have on the beam dynamics. Experimental results (described 
below | indicate that the former approach, described above, gives the closest match to measured responses. 

For purposes of comparison, several experiments were conducted in testing response characteristics of the apparatus. An open 
loop frequency response was found by applying a sinusoidal system input torque (varying the motor current), and recording 
measurements of the tip position; the procedure is similar to that employed in 14). Data was taken over the range 0.2 Hs to 13.0 
Hr. ill steps of 0.1 flz. and the results are shown in Table 3. The system poles and zeros were found by noting the frequencies 
which produced maximum and minimum tip deflection, respectively. An inherent assumption in this technique is that the damping 
of the beam is very small (this fart was experimentally verified in an independent study 19]). The damping ratio calculations 
represented m the table are based on the assumption that excitation near a modal frequency will result in the response showing 
primarily only that particular modal frequency. 


Table 2: FEM Results 

I 

Mode Frequency ( Hertz l Table 3: Frequency Response Data 


1 

2.0091 

Minimum Tip Response Maximum Tip Response 

Damping H. 

- 

8.2509 

1.0 Hi 1.2 Hz 

0.1.19 


23.1187 

KM Hz 7.fi Hz 

0.050 

t 

lti.5ti“7 

It . 1 Hz 12.0 Hz 

0.008 

*» 

79 5211 




The open loop stop response (in position) of the arm was found by rotating the motor shaft through an angle of 10 degrees and 
measuring the tip d-Hertion from its nominal value I initial point). After this maneuver the motor holds the new position (that 
is. is servomg I since the local feedback loop is active. Figure 2 illustrates a plot of the step response. Note that while the torque 
IS applied at the hub at time / - 0, the tip deflection response is delayed by approximately Id milliseconds and. in fact, initially 
moves in the direction opposite that of the hub rotation. The step responses indicates a settling time of about one minute. A 
fast Fourier transform of the data allows ,ear identification of the first two modal frequencies; these o, cur at 1.18 Hz and 7.5 
Hz. respectively. Figure t shows the result of the FFT for the tip position in the step response test. We note that the rigid body 
mode i dr component I due to the pinned joint has been subtracted out of the FFT plot for clarity. 

3. Control Analysis 

/ 1 Pm hi i m Form ul, i/iofi 


Consider again the single link flexible manipulator system described above, redrawn in Figure 1 The displacement .if any point 
along the arm is given by the hub angle (ft/) and the deflection otr,/) measured from the line o>. We assume that only transverse 
vibration is present and that the deflection due to this vibration is small. Let L be the arm length so that in general terms 

ytr.t) = .jfr./l - rtflf I . (11 
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STEP RESP0NSE 


FFT 0F TIP P0SITI0N 



•. i. »». it. m. it. ». m. 

TIME IX 0.009 SECONDS) 
Figure 2 


M l». M. M. «•. M. M. »*. «. 

FREQUENCY (X O.OQS HZ) 


' * Figure 3 

where Ell r) i» (lie elastic stilfnrii. .4 i» the cross-sectional area, and in , is the iiuu density. For the mechanic*! configuration 
under consideration. I 'Jl inuit satisfy the boundary conditioni 


oill.M =0 : £71x1 


- t - 1 '*-* ■ mogi •< • ««*>g! •• • 

(■# t*L * m L 


where T is the turi|iie at the hub and In i* the actuator inertia. Accordingly, (2) may be put into the familiar form for the 
generalized modal coordinates q( I ( as 

A/q - 0q * A'q = Bf . q - i*,«, <j.; r , (4) 

where \l is the mass matrix. A is the stiffness matrix, and D contains terms associated with the damping. For position and 
velocity measurements mi the y direction) the solution to (4) is approximated by 

yir.n ; JslOt'tfx) (5) 

fc •;*» 

for tlw ri^ruUnu turns imlmlm* th«* rigid Ixwly niodr at k - 0. and the J h [l\ are choxeti to minimize the mean -square error 

upon Mih'stitntion into I 1 ». 


t nd**r tin* unitary t ran '.for mat ion th»*ii. w#» um* the notation q now to represent the state in the state variable formulation of 

q 0 / q 0 . f a 1 

q j [ (i : * r D* q * ♦ r fl f ' ' M ' 01 ' ( [q (6) 

where + r /i 'l> i-. .in 'i • ii diagonal matrix containing the squares of the modal frequencies W|.^. along the diagonal. 

.iiul 'I'i i- the niiMsiireiuent vet, ir I pon i (* i ."ran get lie ill . vve may write 1 1> t in the manner 


■/.. 'I 1 

I) I) 


I) 1 


•I- 1 




fs » I ■/„ 0 

.*)'*, I ~' Jr \ * “C.e's . Ifn , , A, 

We note that in ( > I the r.gid body mode lias been included. Since the only control input for this example is the torque, then 
/ T Finally, the matrices (' and + r W are given by 

,, /- 0 fit Z- * 0 (,'slit 0 /■ r ,r 

0 1 0 -*:l» ... 0 ■ + a M l -■'»» - («) 


I. 0 c-ilf.1 
I) I 0 
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i.t Parmmtltr EtlimaUo «i 


The fundamental issue in lh« mathematical formulation of flexible mechanical structures lie* ia the (act that tack distributed 
parameter systems mutt be identified (controlled) with only a limited number of teuton ( actuators), liut-'-d, (or the analysis of 
the tingle-link flexible arm we typically consider hub actuation only, and tip position and/or hub velocity measurements to be 
employed in modeling and feedback control. Moreover, without reliable models for the control design the analysis becomes even 
more difficult. Philosophically, there are several different views to take in the control design. One approach is to construct a 
controller that is reasonably robust in the presence of modeling uncertainties and spillover, and yet simple enough in structure 
to be easily implementable (for example the variable structure control approach [30]). Another approach is to perform system 
identification exercises to model the system as accurately is possible prior to control design. A third approach is a combination 
of tbe first two: estimate the system parameters on-line (in the dosed loop) and base the control design on the resulting model. 
This last viewpoint is often refered to as Self-Tuning Adaptive Control (STAC). 

In tbe STAC approach, the manipulator dynamics are represented by linear discrete-time models, affording tbe primary advantage 
that the controller design is inherently digital in nature. In the application to flexible structures, tuning parameters include 
combinations of the damping and modal frequences, or some combination of other free parameters which make up the manipulator 
model. Our approach to the parameter estimation problem involves recursive least squares methods with covariance resetting. 
That is, in order to maintain a fast overall convergence rate, the covariance of the estimates is reset at regular intervals in the 
algorithm. Such a scheme is particularly attractive for the manipulator control problem due to the time-varying nature of the 
tuning parameters during slew maneuvers and varying payload exercises. Experimental studies of the parameter estimation and 
STAC approach for the arm described above are presently underway. In the following we present simulation results which indicate 
avenues to pursue regarding implementation. 



Figure 4: Deflection Slew angle tf(<) 
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Figure 5: Parameter Estimation Simulation Results 
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The dementi of the C matrix can be found according to jH| 


*<M£<orfel (9) 

and the (dVi/df )(0) are solved from a system of nonlinear equations. The modal frequencies can be computed a priori, or identified 
as discussed in the previous section. We model the system as a stochastic ARMA process and excite a fourth-order model of the 
arm with a white noise input. Such a representation allows a delay (in tip position response, as observed in the actual system) 
to be inserted into the model. Using a aero-order hold circuit in the model and sampling the simulated FEM model, the AR and 
MA parameters converge to their nominal values as depicted in the sample plots of Figure S, which shows time histories of one 
AR parameter and one MA parameter. Values for damping coefficients are then calculated from these parameters. These results 
are not useful in closed-loop control however, due to the length of time for convergence to the true parameters. Note also that a 
primary difficulty results because of the approximate pole-sero cancellation in the system model (indicated by the spike at about 
0.7 se.onds). Slightly better results are obtained if the rigid mode is removed from the model, which corresponds to exciting the 
un/onrof system with an initial disturbance. Such an exercise is possible since the motor inertia is considerably greater than that 
of the arm. 

Prior to actual experimentation on the arm, several modifications must be investigated. For example, simulation studies for this 
and t-.lier example systems have indicated improvement for different resetting intervals; for details, the reader may wish to consult 
[21 j. Also, simple digital low-pass filtering of the measured variables has produced improved performance of the parameter 
estimator. For control purposes the simulations have shown that an algorithm which turns on the control after allowing the 
estimator to run for a short period of time (for example, as illustrated in the simulations, about 1 to 1.5 seconds) will achieve 
the control objective. However, we are presently pursuing ways of improving the time to convergence in the closed loop with 
approaches using state feedback. 

J.J Output Firilhnri unit Frequency Shaping 


Generally speaking, high dimensionality and multiplicity of inputs in large-scale systems such as flexible mechanical structures 
leads to complex centralized controller schemes. One solution to this problem is to simplify the structure of the model via decom- 
position into subsystems with associated subcontrollers in a decentralised output feedback formulation. Moreover, centralized or 
decentralized output feedback is one of the more straightforward algorithms, from the viewpoint of implementation, for the control 
of flexible mechanical slructur*s: see, for example, [22,2.1.201. 

For the problem of single-link flexible manipulator control, where only hub actuation is employed in the control action, the output 
feedback approach to controller implementation is centralized in nature. The problem of spillover is, however, a critical issue to 
consider in the design. In order to minimize the effects of spillover, we consider a frequency-shaped cost functional [24) , where 
penalties aie assigned to the truncated modes and high penalties are assigned to the high harmonics at the input in order to 
minimize the effects due to excitation of the residual modes. 

We consider the cost functional J to be minimized as formulated in the frequency domain utilizing Parseval’s Theorem. With 
infinite time horizon, such a cost is written in the manner 


J = r .Y*(;w)<?0ui).Y0'w) + r*(;w)R(;ui)T(;ui)|du- , 

J - TO 


where Y(j^| corresponds to the system state (as in (fi)) and T(ju>) the input (torque) of the system. For implementation of such 
a scheme, consider the diagram of Figure 6, where the parameters A’|, A'j, A's are solved for in the minimization of (10), and the 
filter pole location I ■> I is dependent on the system dynamics. In the example under consideration, Q(j u>) is the system matrix and 
Rlj*t I - lyw -■> /( - - i ). 

I’nder this formulation, the open -loop state variable representation of the system has the form 
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Figure 6: Output Feedback Scheme 
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where the new measurement v(f) now includes the torque T. Incorporating the system of (11 )— (12) into the output feedback 
structure of Figure 6 (and subsequent solution of the corresponding Lyapunov equation) allows the off-line calculation of the 
feedback gains from minimization of (10) by an appropriate nonlinear optimisation routine. 


We consider now a simulation of the flexible arm system, using a five-mode model from the FEM as the “truth model”, from which 
measurements are taken and fed back in the output feedback scheme. The controller design for this example is based on the reduced 
system of rigid mode plus first flexible mode, and the resulting control is then tested against the full-order truth model to illustrate 
the effects of the frequency weighting approach in reducing spillover. A conjugate gradient method is used in the optimization 
portion of the design, and the final values obtained for the control law (with 7 = 4) are A'i = - 1 10.09, A'i = - 1 1 1.53, A'i = -88.83 
(feasible values for the system under consideration) for the cost which reached a minimum after approximately 3000 iterations. 
The results using this controller are illustrated in Figure 7 for a step input torque; this applied input is such that the tip rotates 
through a small angle (of less than 5’, in terms of the rigid position). The values for torque begin at zero, that is, the dc component 
is subtracted out. In the simulation, the response settles in about ten seconds, whereas the free response decays after about one 
minute due to damping included in the model. 


4. Structure-Mounted Proof-Mass Actuation 


Sirce the large-angle slewing problem is complicated due to the flexibility effects inherent in the structure to be slewed, one is 
naturally led to investigate the possibility of relegating, at least partially, the task of vibration damping to a separate sensor- 
actuator pair and associated feedback loop. To this end, we have been investigating utilization of a structure-mounted momentum- 
exchange device mounted near the tip of the single-link manipulator. The practicality of such active vibration damping in a 
robotic environment is closely coupled to the availability of lightweight and effective devices. The device we have considered in 
our preliminary studies is a proof-mass actuator developed in our labs to study active vibration in space based flexible mechanical 

structures. 

Non-ground referenced linear actuators are not yet widely available on the market, and this fact led to an in-house development; 
a general view of the device as mounted on the arm is shown in Figure 8. The device is built around a linear motor manufactured 
by the Kinico division of BEI Motion Systems which has a total mass of 25 grams, and can deliver a peak force of 2.2 N. The 
coil (solenoid) is rigidly mounted to a beam clamp which fixes the actuator to the arm. Also connected to the clamp is a rigid 
aluminum bracket which supports the springs. The proof mass is coupled to the framework through the springs, which are in 
turn coupled to the framework with adjusting screws so that their rest tension and proof mass rest position can be controlled. 
There is sufficient adjustment so that springs of different length and stiffness constant can be accomodated. The springs provide 
a restoring force for the proof mass and transfer force to the structure. A hanger was also mounted to provide strain relief for the 
feed wires. 

The proof mass consists of a rectangular steel ring with a central steel member. This central member passes through the coil and 
restricts motion to a single axis. Samarium cobalt magnets are fixed to the top and bottom inside edges of the ring (adjacent to 
the coil) so that the interaction of the permanent magnetic field with a current in the coil results in a force on the proof mass. 

Details ot the development of the dynamic model of the actuator may be found in (17|. The net force applied to the tip of the 
arm (where the actuator is mounted) may be given as 

/ =2kz - Kri. + B- , (13) 
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Figure 7: Output Feedback Example 



mi + B: + 2kz x Kfi + mj , (14) 

where m is the proof mass, y the displacement of the structure at the point of actuator attachment, / the force acting at that 
point, : the relative displacement of the proof mass, B the viscous damping coefficient, hr the motor force constant, and i, i« are 
the input and armature currents, respectively. The actuator constants taken from the data sheets which accompany the individual 
components, are m = 0.0207 kg, fc = 262.7 N-m' 1 , A> = 1.112 N-ampere-*. The incorporation of the above actuator creates a 
second feedback loop to which the task of vibration damping is relegated. The two control loops (for slewing and for vibration 
damping) can be both designed and implemented in a decentralized manner. Note that y includes the displacement due to both 
the rigid body mode and the flexibility (see (l)). The principle of relegation implies that we design the feedback control only for 
the latter portion. To this end consider a vibration damping loop for only the first mode, such that the relevant expression is 

it + = btf i (I*) 

where -/| is the natural frequency of the first mode and b, is an influence factor determined from the mode shape at the point of 
interest. Acceleration feedback can be used from the co-iocated accelerometer and a simple PI controller has been designed. It 
is evident, however, that the STAC approach or the frequency weighted control approach outlined earlier, can also be used here. 
The incorporation of the more sophisticated design approaches resulting in more complicated controllers will aid in handling more 
than the first vibrational mode. Studies along this direction are presently continuing. 

S. Conclusion 

In this workshop presentation we have described work in progress on modeling, parameter estimation, and control studies for 
an experimental, one-meter single-link flexible manipulator arm. Models have been developed for the apparatus based on finite 
element analysis and experimental verification. These, with the closed-loop parameter estimation procedures described here, and 
subsequent STAC approach for control, are being evaluated on the laboratory arm. 

Under investigation is experimentation involving local proof-mass actuation for vibration control at the tip of the arm, using a 
device developed in the Control Research Laboratory at Ohio State for flexible structures control work. The output feedback 
frequency shaping approach described here may be easily extended to this application, where the formulation is decentralized in 
nature: results on this technique for general flexible structure vibration control will appear in [25j. Finally, various other centralized 
(for tin" case of \nl> actuation only) ami decentralized approaches are currently being evaluated in the laboratory. 
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